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ABSTRACT 

We report the discovery of a minor planet (2006 SQ372) on an orbit with a perihelion of 24 AU and 
a semimajor axis of 796 AU. Dynamical simulations show that this is a transient orbit and is unstable 
on a timescale of ^200 Myrs. Falling near the upper semimajor axis range of the scattered disk and 
the lower semimajor axis range of the Oort Cloud, previous membership in either class is possible. 
By modeling the production of similar orbits from the Oort Cloud as well as from the scattered disk, 
we find that the Oort Cloud produces 16 times as many objects on SQ372-like orbits as the scattered 
disk. Given this result, we believe this to be the most distant long-period comet ever discovered. 
Furthermore, our simulation results also indicate that 2000 OO67 has had a similar dynamical history. 
Unaffected by the "Jupiter-Saturn Barrier," these two objects are most likely long-period comets from 
the inner Oort Cloud. 

Subject headings: Oort Cloud, long-period comets 



1. INTRODUCTION 

In the past three decades, many minor planets have 
been observed on planet-crossing orbits in the outer 
solar system. These objects are transient in nature 
(jTiscareno fc Malhotral 12003 iHolman fc Wisdornl Il993f ) 
and, therefore, must be continuously resupplied from 
more stable reservoirs. Previous numerical modeling 
has shown that both the scattered disk and Oort Cloud 
contain large populations of bodies diffusing into the 
planetary region on Gyr timescales at a rate that can 
account for the planet-cr o ssing populations observed 
today (jLevison fc Duncanl 119971 : IWiegert fc Tremaind 
1199a lEmel'Yanenko et al. 2007). In the scattered disk, 
trans-Neptunian objects (TNOs) orbit the Sun on eccen- 
tric orbits with perihelia taking them to within 10 AU 
of Neptune or less. Driven by the effects of Neptunian 
perturbati ons, these orbits underg o a chaotic evolution 
over time (jDuncan fc Leviso"nlll997| ). Due to this process, 
some of these bodies eventually begin to scatter more 
strongly off of Neptune, and because complete ejection 
by Neptunian encounters is quite unlikely, t he perihelia 
of thes e orbits tend to evolve further inward (iFernanded 
I1980H) . 
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While the scattered disk is typically assumed to be the 
source o f most observed plan et-crossers in the outer Solar 
System ([Duncan et al.ll2004T ). there is reason to expect a 
substantial Oort Clou d contribution as well , espec ially at 
high semimajor axes (|Emel'vanenko et alj [20051) . Oort 
Cloud orbits are continually driven into the planetary 
region via perihelion shifts due to the extern al perturba- 
tions of the Galactic tide and passing star s ([Oortl[l950l : 
lFernandezlll980al ; lHeisler fc Tremaindll986f h When Oort 
Cloud bodies have their perihelia pushed back into the 
planetary region, their orbital ev olution is l argely de- 
termined by their semimajor axis (|HillsllT98l[ ). For low 
semimajor axis orbits (a < 20,000 AU), perihelion migra- 
tion through the planetary region occurs very slowly over 
many orbital periods because the external perturbations 
driving this migration are small compared to the binding 
force of the Sun. Thus, as these objects slowly drift Sun- 
ward, they will inevitably encounter Jupiter and Saturn, 
receiving a gravitational energy kick that ejects these 
weakly bound bodies to interstellar space. On the other 
hand, the perihelia of high-scmimajor axis Oort Cloud 
objects (a > 20,000 AU) evolve much more rapidly, and 
these bodies' perihelia can move from outside 15 AU to 
inside 5 AU in less than one orbital period. As a result, 
before being ejected by Jupiter or Saturn, these objects 
can pass through the terrestrial planet region where we 
subsequently detect them as long-period comets (LPCs). 

Thus, the Oort Cloud can be divided into an outer 
Oort Cloud (a > 20,000 AU), from which all known 
LPCs originate, and an inner Oort Cloud (a < 20,000 
AU) that is not sampled by known LPCs but is theo- 
rized to exist due to mode ls of the comet cloud's for- 
mation (jDuncan et alJ ll987). Because of this, obtaining 
a sample of comets that probes the entire Oort Cloud 
requires detecting LPCs outside the gravitational bar- 
riers of Saturn and Jupiter (q > 15 AU). Due to the 
steep decrease in brightness with heliocentric distance, 
however, detection of these objects is very challeng- 
ing. Until this work only one LPC, Comet C/2003 A2 
(q = 11.4 AU), has been observed with a perihelion be- 
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yond Saturn's orbit (jGleason et all 120031 : iGreenl |2003h . 
Detecting LPCs bey ond ~15 AU w o uld te st Port Cloud 
form a tion models (iDuncan et all 119871 iDones et al.l 
[200l lEmel'Yanenko et al.l 120071 : iKaib fc Quinnl I2008D 
as well as constrain the birthplace environment of 
the Sun dFemanded 11997k iFernandez fe Bruninil 120001 : 
iBrasser et al.l 120061: IKaib fc Quinnl I2008D . For our pur- 
poses, we define an LPC as any Oort Cloud body that has 
been reinjected into the planetary region. These objects 
may or may not show cometary activity. 

Here we report a likely LPC well outside the Jupiter- 
Saturn barrier with q = 24.2 AU. Section 2 discusses 
our observations of this object. Section 3 is devoted to 
modeling the production of similar objects from both the 
Oort Cloud and the trans-Neptunian scattered disk. In 
Section 4, we use our simulations to make the case that 
2006 SQ372 is indeed from the inner Oort Cloud rather 
than the scattered disk. 

2. OBSERVATIONS 

2006 SQ372 was initially discovered in images taken by 
the SDSS-II Supernova Survey (Frieman et al. 2008; see 
also Fukugita et al. 1996, Gunn et al. 1998, Gunn et al. 

2006, York et al. 2000). This survey repeatedly scanned 
a 300 sq. deg. region (SDSS Stripe 82) centered on the 
celestial equator in the southern Galactic hemisphere for 
three months out of the year (Sept - Nov) from 2005- 

2007. This area also passes through the Solar System 
ecliptic plane, spanning ecliptic latitude —20 deg < [3 < 
20 deg and ecliptic longitude —58 deg < A < 58 deg. A 
search through these data by Becker et al. (In prepara- 
tion) attempted to link all single-epoch transients into 
orbits consistent with Kepleria n motion around the Sun, 
using the linking softw are of iKubica et al.l (120071) an d 
orbit fitting software of iBernstein fe Khushalani ( 2000). 
Nearly 50 trans-Neptunian objects were discovered in 
this effort, including 2006 SQ372. 

2006 SQ372 was first discovered as a single-epoch 
transient on the night of Sept 27, 2006 at Right As- 
cension 21:09:20.3 and Declination +01:09:36 (J2000), 
at an apparent magnitud e of r = 21 .5 9 ± 0.06 in 
the magnitude system of iLupton et al.l (|1999fl . The 
data are ast rometrically calibr ated using the methods 
described in iPier et ahl (|2003l ) and have been photo- 
metrically r e calibr ated using the methods described in 
llvezic et ahl (|2007l ). Subsequent observations through- 
out the 2006 observing season recovered this object 
on 9 additional dates as it moved through Stripe 82. 
These observations were used to calculate an orbital 
ephemeris, and to precover additional data from 10 
nights in 2005, as well as 7 subsequent nights in 
2007. The orbital ephemeris is available through the 
minor planet data center, and detailed observational 
data can be provided on request and is also posted 
at !htt p://www.boulder.swri.edu~buie/kbo/astrom/| 
06SQ372.html. 

The current orbital arc now spans 2.9 years. Using 
the Jet Propulsion Laboratory (JPL) Horizons soft- 
ware 1 , the barycentric orbital elements for this system 
at epoch 2008 May 14.0 are semi-major axis a = 796 ± 3 
AU, eccentricity e = 0.9696 ± 0.0001, and inclination 
i = 19.46471 ± 0.00008 degrees. For comparison, the 

J http:/ /ssd.jpl.nasa.gov/?horizons 



Minor Planet Center (MPC) software provides a best- 
fit orbit of a/e/i = 1057/0.9771/19.46052 evaluated in 
heliocentric coordinates. For objects at such high ec- 
centricity, barycentric coordinates are more stable than 
heliocentric coordinates, and orbital elements are evalu- 
ated with respect to the Solar System barycenter for the 
remainder of this work. 2006 SQ372 reached its perihe- 
lion of 24.1667 ± 0.0004 AU on 2006 Aug 12, and has an 
orbital period of 22,466 years. 

In Figure 1, we compare the perihelion and semimajor 
axis of 2006 SQ372 to all other TNOs and Centaurs listed 
in the MPC. We see that the semimajor axis of this object 
is larger than any other known body with perihelion in 
the outer Solar System. The next two largest semimajor 
axes are those of 2000 OO67 and Sedna. While Sedna's 
large periheli on makes it dynamically iso lated from the 
outer planets (|Morbidelli fc Levisonll200l . 2000 00 67 is 
in a Neptune-crossing orbit, making it most similar to 
2OO6SQ372. 
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Fig. 1. — Plot of q vs. a for all listed TNOs and Centaurs in the 
Minor Planet Center. The orbits of Sedna, 2000 6 7, and 2006 
SQ372 are shown with diamond, triangular, and square data points 
respectively. All other observed objects are displayed with smaller 



3. NUMERICAL METHODS 

All numerical work in this paper is done using the 
SWIF T RMVS3 integr ation package (iLeyison fc Duncan! 
119941) modified as in IKaib fc Quinnl (|2008l ). and in- 
cludes the gra vitational effects of th e Sun, giant planets, 
assing stars (Rickman et al.ll2008D . and Galactic tide 
Levison et al.l 20011 ). We perform three different simu- 
lations to investigate the dynamics of 2006 SQ372. The 
first is a simulation of 10 3 test particles cloned from the 
observed orbit of 2006 SQ372 to determine how this orbit 
will evolve in the future. This first simulation motivates 
our other numerical work, which we use to constrain the 
dynamical origin of this object. To determine the most 
likely origin scenario for 2006 SQ372 , we run two different 
simulations: one modeling the production of similar bod- 
ies from the scattered disk and the other from the Oort 
Cloud. It must be noted that this approach assumes that 
the scattered disk and the Oort Cloud are the only two 
populations capable of producing planet-crossing orbits 
in the outer solar system. While previous work has shown 
that classical KBOs and resonant TNOs can also evolve 
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on to Neptune-c rossing orbits (jDuncan fc Leviso"nlll997l : 
iMorbidellilfTOQl . the occurence of this i s substantially 
less f requent than in the scattered disk ([Duncan et al.l 
2004). Our three different simulations are described in 
the following subsections. 

3.1. Preliminary Clone Simulations 

With a perihelion inside Neptune's orbit, the orbit of 
2006 SQ372 appears to be transient in nature. To verify 
this we first integrate 10 3 clones of this object. These 
clone orbits are generated from a multivariate distribu- 
tion using the covariances between orbital elements. We 
find that the population of clones in planet-crossing or- 
bits has a half-life of -180 Myrs. After 1 Gyr, only 2% 
of our clones are still found in Neptune-crossing orbits, 
while 62% have been ejected from the solar system (r > 
200,000 AU), and 36% reside in the Oort Cloud (q > 45 
AU). 

Because 2006 SQ372 is on a planet-crossing orbit not 
stable over the lifetime of the Solar System, it must 
have previously belonged to another class of objects. We 
believe the Oort Cloud and scattered disk are the two 
most likely source populations for 2006 SQ372, as both 
of these efficiently pla c e bod i es on planet-crossin g orbit s 
(|Wiegert fc Tremaind 119991: iLevison fc Duncan! Il997f ). 
An Oort Cloud origin would imply that torques from 
passing stars and the Galactic tide reinjected its perihe- 
lion into the planetary region after which planetary en- 
counters decreased its semimajor axis. However, it is also 
possible that 2006 SQ372 evolved from the scattered disk. 
In this scenario, the perihelion of 2006 SQ372 chaotically 
diffused inward before its semimajor axis was inflated via 
planetary encounters. 

3.2. Scattered Disk Production 

To determine the likelihood that SQ372 is from the 
scattered disk, we simulate the production of similar or- 
bits fro m this population. To do this, we replicate the 
work of ldi Sisto fc Brun ini (2007|) who modeled the pro- 
duction of Centaurs from orbits cloned from observed 
scattered disk objects. Particles are removed from our 
simulation if they reach r > 6,000 AU or collide with a 
giant planet or the Sun. 

We begin our simulation with 2,500 test particles and 
the four giant planets on their present orbits and inte- 
grate for 4.5 Gyrs. The initial conditi ons we use to assign 
test p articl e orbits are very s imilar to Idi Sisto fc Brunini 
(2007) and lFernandez et all (|2004l ) and are based on the 
orbital distribution of observed SDOs. To first obtain a 
sample of real SDO orbits, we classify any observed TNO 
with 30 AU < q < 40 AU and a > 50 AU as an SDO. 
Additionally, we also include observed TNOs with 30 AU 
< q < 40 AU, 40 AU < a < 50 AU, and e > 0.2. This 
selected set of real orbits then forms the basis for our ini- 
tial perihelion distribution, which we generate as follows. 
First, we weight the perihelion distribution of real orbits 
by the fraction of orbital period each object spends in- 
side 45 AU, an approximate zone of observability. Using 
this weighted distribution, we then assign perihelia to our 
test particles. Next we assign semimajor axes between 40 
AU and 200 AU to particles by randomly drawing from a 
distribution proportion al to a~ 2 . We choose this power- 
law distribution because fdi Sisto fc Brunini (|2007f ) found 
it to be a good fit to their debiased distribution of real 



SDO semimajor axes. To assign orbital inclinations, we 
us e the o b serve d SDO inclination distribution obtained 
in iBrownl (12001 of 

F (i) oc sin i exp -i 2 /2of (1) 

where we choose o~i — 12° as in Idi Sisto fc Bruninil 
(|2007| ). Lastly, arguments of perihelion, longitudes of 
ascending node, and mean anomalies are assigned ran- 
domly. 

Because we are interested in the distribution of orbits 
penetrating the planetary region from the scattered disk, 
we clone a particle 20 times if it evolves to an orbit with 
q < 30 AU during the course of our simulation. This 
cloning is done by shifting each cartesian coordinate by 
a random number between ±1 x 10 -6 AU. These clones 
are then evolved until collision with the Sun or a giant 
planet, or ejection (r > 6,000 AU). 

We choose to remove particles with r > 6,000 AU be- 
cause we do not want to include any SDOs that move to 
the Oort Cloud and then subsequently evolve to planet- 
crossing orbits from there. However, doing this cuts off a 
possible dynamical pathway where a temporarily random 
walks above 3,000 AU and is subsquently pulled back 
down to SQ372 values while q never leaves the planetary 
region. To evaluate the importance of this pathway, we 
continue to integrate a subset of the particles removed at 
r = 6,000 AU. Based on these integrations, we find that 
SQ372-like orbits are attained about an order of magni- 
tude less frequently than for particles yet to reach r > 
6,000 AU. Thus, our simulation does capture the dom- 
inant mechanism producing 2006 SQ372 analogues from 
the scattered disk. 

3.3. Oort Cloud Production 

To evaluate the rate the Oort Cloud produces orbits 
similar to 2006 SQ372, we dynamically evolve 10 6 Oort 
Cloud particles for 1.4 Gyrs under the influence of the 
giant planets, Galactic tide, and passing field stars. To 
simulate this large number of particles efficiently, we 
ha ve used the var i able t imestepping procedure outlined 
in iKaib fc Quinnl (|2008h with an extra level of larger 
timesteps. The base timesteps are t =0.25, 10, and 50 
yrs for r < 300 AU, 300 AU < r < 400 AU, and r > 400 
AU respectively. 

To choose our initial particle orbits, we use orbital el- 
ement distributions from the ~ 10 3 -particle Oort Cloud 
that e xists at the end of the 4.5 Gyr large control simula- 
tion of lKaib fc Quinnl (|2008h . First, we assign semimajor 
axes by randomly drawing from the a cumulative distri- 
bution function (CDF) from our formation simulation. 
In this formation simulation a ranges from 911 AU to 
~10 AU. Next we divide our formation simulation into 
five separate ranges of a, each containing an equal num- 
ber of particles. For each a-range, we then construct 
CDFs for J/L, J z /J, fi<3, and u>g, where L is the first 
Delaunay momentum (wGMqo), J is the total angular 
momentum, J z is the z-component angular momentum 
in the galactic frame, VLg is the longitude of ascending 
node in the galactic frame, and is the argument of 
perihelion in the galactic frame. These CDFs are then 
used to randomly assign e, i, fl, and u> for the new ini- 
tial orbits in each a-range. Finally, mean anomalies are 
assigned randomly from a uniform distribution. 
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The goal of this procedure is to build a 10 6 -particle 
Oort Cloud that looks very simi l ar to the ~ 10 3 -particle 
cloud formed in iKaib fe Quinnl (|2008f ) . This is of par- 
ticular importance in the inner Oort Cloud where or- 
bits are less isotropized. Due to the very long Kozai 
cycle times for the galactic tide inside a < 5,000 AU, 
the extreme inner Oort Cloud should contain signifi- 
cantl y more orbits with 17 > t han with q < ((Kozail 
Il962t iHeisler fc Tremaing |l986). Failing to replicate 
this effect would enhance the rate that extreme in- 
ner Oort Cloud objects enter the planetary region and 
could artficially in crease the production of o rbits simi- 
lar to 2006 SQ 372 . IHeisler fc Tremaind (|1986h show that 
q oc —y/qa 2 sin 2ojq. Because we apply a different ujq dis- 
tribution for each range of a, our 10 6 -particle simulation 
does capture this bias at small a. Moreover, it is likely 
that much of the inner Oort Cloud was randomized by 
perturbations of an early star clust er environment not 
included in the lKaib fc Quinnl (|2008f ) cont rol simulation, 
and t his would greatly diminish this q bias (jBrasser et al.l 
2008). For this last reason, we can assume that the ini- 
tial conditions for our simulation are conservative in this 
respect. 

Because planet-crossing orbits typically result in ejec- 
tion from the Solar System, there should be a constant 
deficiency of these orbits within the Oor t Cloud. How- 
ever, the orbital distributions from the Ka ib fc Quinnl 
(2008) Oort Cloud formation simulation are too coarse 
to resolve this deficiency. As a result, this part of the 
Solar System is initially overpopulated with Oort Cloud 
orbits. To correct for this effect and allow the planets to 
eject these bodies, we evolve our comets under the grav- 
itational influence of the Sun, four giant planets, passing 
stars, and the Milky Way tide for ~1 Gyr. In the fi- 
nal 400 Myrs of our simulation, external perturbations 
rather than our initial conditions control the population 
of planet-crossing orbits at all Oort Cloud semimajor 
axes. 

4. SIMULATION RESULTS AND DISCUSSION 
4.1. Comparison with Previous Work 

To confirm that our simulations model the dynamics of 
the scattered disk and Oort Cloud properly, we compare 
our generated dynamical populations with past work. 
First, in the case of the scattered disk populatio n, we 
replicate the analysis of Idi Sisto fc Bruninil (|2007f ) who 
model Centaur production from SDOs and obtain a time- 
weighted perihelion distribution for our Centaur popula- 
tion. Upon comparing the two perihelion datasets, we 
find that a K-S test gives a 97% probability that they 
were drawn from the same distribution. Furthermore, 
we can also examine the mean lifetime of Centaurs in 
our simulation. When we apply th e q < 5.2 AU and 
a > 1 000 AU removal rules used in Idi Sisto fc Bruninil 
(|2007f ) to our data, we obtain a mean Centaur lifetime of 
77 Myrs, which is comparable to their value of 72 Myrs. 

In the case of our O ort Cloud model, we c ompa re our 
simulation to that of lEmeP Yancnk oet al.l (|2007l ) who 
modeled the production of comets from the Oort Cloud. 
By measuring the rate that dynamically new comets en- 
ter the inner 5 AU of the solar syste m in their simulations 
and as suming a real rate of 10 vr~ 1 . lEmel'Yanenko et all 
(|2007f) calculate that there must be 4.5-5.5 x 10 11 bodies 



in the Oort Cloud greater than 1 km in diameter. We 
perform an analogous calculation with our Oort Cloud 
simulations and determine that the Oort Cloud must 
contain ~5 x 10 11 bodies, which is in excellent agree- 
m ent. Even more enc o uragi ng is that the simulations 
of lEmel'Yanenko et al.l (|2007h predict objects similar to 
2006 SQ372 passing through the outer solar system from 
the Oort Cloud. 

4.2. 2006 SQ 372 Production Rates 

The goal of our numerical work is to determine the rel- 
ative rates that 2006 SQ372 analogues are produced from 
the Oort Cloud and scattered disk. To do this, we sample 
the orbital elements of all particles at discrete timesteps. 
For the scattered disk simulation, we sample our cloned 
orbits every 10 5 years, and our Oort Cloud simulation is 
sampled every 10 6 years. To minimize the effects of our 
initial conditions, we only use the latter portions of each 
simulation. For the scattered disk, we analyze the last 
3.5 Gyrs, whereas in the Oort Cloud simulation we only 
use the final 400 Myrs. We also choose these time lengths 
because they yield roughly equal numbers of orbits with 
q < 30 AU. Using this time-sampled data, we then gener- 
ate a two-dimensional histogram of orbits in g-log a space 
for both datasets. The grid spacing we use is 2 AU for 
q and 0.278 for log a. Tracing the particles' evolution as 
they traverse through our q-log a histogram then allows 
us to construct a residence time distribution for orbits 
leaking into the planetary region from each reservoir. 

Once we have compiled orbital residence time his- 
tograms for both of our simulations we must now find 
a way to normalize each histogram bin population to 
match our own solar system. In the case of our scattered 
disk simulation, we can use the number of Jupiter Fam- 
ily Comets (JFCs) produced to normalize our histogram, 
whereas for our Oort Cloud simulation we can normalize 
our dataset to match the rate new LPCs are produced. 
Because the observed JFC population is dominated by 
bodie s greater than 1 km in diameter (|Tancredi et al.l 
2006), we impose a similar minimum diameter (D > 1 
km, or Hio < 11) on the observed flux of LPCs used 
in our normalization routine. Unfortunately, numerous 
parameters used to generate these normalization factors 
are uncertain by up to an order of magnitude. In the 
following, we explore a range of normalization factors for 
each simulation. 

We begin by discussing normalization our scattered 
disk simulation. In our own Solar System, there are 
250 known active JFCs (2 < Tj < 3, q < 2.5 AU). 
Given their proximity to Earth and their bright abso- 
lute magnitudes, we can take this as an observationally 
complete sample. In addition, however, there should be 
a second population of dormant J FCs. Based on mod- 
els of physical aging, Uewittj (|2004h estimates that there 
are about twice as many dormant JFCs as active ones. 
In contrast, by modeling the diffus i on of JFC inclina- 
tions with time, iLevison fc Duncan] (|1997f ) argue for a 
3.5:1 dormant to active JFC ratio and possib ly as high 
as 7:1. Lastly. iFernandez fc Morbidellil (|2006h find a ra- 
tio between 1 and 3 by studying dormant JFCs found in 
NEO observations. I Yolk fc Malhotral (|2008h . who derive 
the size of the scattered disk population based on known 
JFCs, adopt a ratio of 2 because of the agreement be - 
tween IFernandez fc Morbide lli ( 2006) and lJewittl (|2004h . 
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For the same reason, we choose this value as the most 
probable one as well. This implies a total number of 
750 JFCs in the current solar system. Because we find 
546 instances of JFC orbits in our time-sampled data, we 
multiply the population in each of our histogram bins for 
the scattered disk simulation by a factor of 750/546. 

In the case of our Oort Cloud orbital residence time 
histogram, we use the production rate of km-sized LPCs 
to determine a normalization factor. Unfortunately, the 
observed flux of new LPCs in t he inn er solar system is 
not well-constrained. lEverhartl (|1967al lbl) give a value of 
~10 dynamically new comets per year inside 5 AU, and 
this is used in many nu merical works t hat model the Por t 
Cloud popul ation e. g.. iHeislerl (Il990l); IWeissmanl (|1996f ); 
iDones et all (|2004l) : lEmel'Yanenko et al.1 (|2007n . How- 
ever, using only LPCs observed in the LINEAR survey 
(wher e observationa l biases are much better character- 
ized) , Jfcincii |2005|)finds a shallower size distribution 
than lEverhartl ( 1967bf ) and argues for a flux of 4 new 
LPCs/yr inst ead. Furt h er an alyzing the LINEAR LPC 
observations, iNeslu san (2003) finds an even lower new 
LPC flux due to different discover y efficienc i es bet ween 
new and returning LPCs. Although lNesrusanl (|2007f ) does 
not directly spec ify an LPC flux, the paper states that 
the lHeisierf (|1990f) assumption of 2.1 new LPCs/yr inside 
1 AU is likely an order of magnitude too large. In our 
simulations, this would correspond to an observed flux 
estimate of 1.5 new LPCs/yr inside 5 AU. To be con- 
servative, we adopt this last flux estimate. In our Oort 
Cloud simulation, we find that 9 dynamically new LPCs 
pass through the inner 5 AU of the solar system every 
Myr on average. As a result, we multiply all of our Oort 
Cloud histogram bins by [1.5/(9 x 10" 6 )]/400. The fac- 
tor of 400 is included because we have performed 400 
different time samplings of orbital elements, and we are 
normalizing to a rate of comet production rather than a 
steady state population as in the case of JFCs. 

In addition to the uncertainties in normalizing our sim- 
ulations to JFCs and LPCs, our results are also affected 
by the uncertainty in the inner Oort Cloud population 
size. The reason the inner Oort Cloud population size is 
poorly known is because no LPCs near Earth come from 
this reservoir. Because of this, its uncertainty has no 
impact on our overall dataset normalization, but it may 
have significant consequences on SQ372 production rates 
since inner Oort Cloud bodies with q > 15 AU are not de- 
flected by the Jupit er-Saturn barrie r . The simulation we 
use is based on the iKaib fe Quinnl {2008) control simu- 
lation, which predicts a 1.5: 1 inner-to-ou t er Po rt Cloud 
population ratio. Although IDones et al.l (|2004f) finds a 
slightly lower ratio of 1:1, most simulations modeling 
Oort Cloud for mation predict higher ratios, with some 
as high as 10:1 (jBrasser et alj|2006t lEmel'Yanenko et al.l 
l2007tlKaib fc Quinnl 120081 : IBrasser et alT 2008) . Because 
the inner Oort Cloud population was likely enriched by 
strong pe rturbations of the S un's birthplace environ- 
ment, the IKaib fc Quinnl (|2008l ) control simulation prob- 
ably underestimates the inn er Oort Cloud popul ation 
size. Based on the result s of IKaib fc Quinnl (|2008j ) and 
IBrasser et all ()2006l . l2008l) . we believe that a 3:1 inner-to- 
outer Oort Cloud population is reasonable. To account 
for the full range of possible population ratios, we weight 
particles with an initial a < 20,000 AU by a given factor 
to model more anemic or enriched inner Oort Clouds. 



For instance, we weight these particles by a factor of 2 
to model the orbital distribution of our chosen inner-to- 
outer Oort Cloud ratio of 3:1. 

With both of our histogram datasets normalized to the 
current solar system and an inner Oort Cloud model cho- 
sen, we now have orbit distributions through the plane- 
tary region for bodies originating from both the scattered 
disk and Oort Cloud. Next we divide the Oort Cloud 
histogram dataset by the scattered disk one. The result 
gives the ratio of orbits produced from the Oort Cloud 
to the number produced from the scattered disk in each 
histogram bin. Now for an orbit with a given q and a we 
have a means to estimate the probability that it is from 
the Oort Cloud vs. the scattered disk. 

Before proceeding further, it is necessary to determine 
how accurate this method of mapping orbit production 
is. For instance, if a particle is improbably scattered into 
a relatively stable orbit, its orbit will be sampled over and 
over for the entire simulation. As a result, our map will 
be artificially distorted to reflect a high production rate 
for that particular orbit, even though the chance of the 
original scattering may be quite small. To guard against 
this effect, we attempt to simulate a large enough num- 
bers of particles so that this type of error is diminished 
to an acceptable level. To test whether we have achieved 
this, we subdivide our particles into 4 groups, treat each 
one as a separate simulation, and then compare the re- 
sults. We find that the different probability ratios for 
the orbital bins near 2006 SQ372 have a standard devi- 
ation of 16%. We can take this standard deviation as a 
conservative estimate of our real error. 
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Fig. 2. — Contour plot of the ratio of Oort Cloud origin probabil- 
ity to scattered disk origin probability. When generating this plot 
we only consider orbits with inclinations between 10° and 30° . The 
'X' marks the orbit of 2006 SQ372, and the triangle marks the or- 
bit of 2000 OOe7 The dashed contour lines mark origin probability 
ratios of 10:1 and 1:1. 



Our resulting probability ratio map is shown in Figure 
2. For this analysis, we use the normalization factors and 
inner Oort Cloud population that we have determined 
above. In this figure, only orbital inclinations between 
10° and 30° are included, as these most closely corre- 
spond to the true inclination of 2006 SQ372. It should 
be noted that the scattered disk contribution decreases 
with higher inclinations since it is concentrated toward 
the ecliptic much more than the Oort Cloud. It can be 
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TABLE 1 

Parameters used in generating the Oort Cloud to scattered disk probability ratio and the resulting ratios. Columns in 
order are: (1) the number of active jfcs in the current solar system, (2) the ratio of dormant to active jfcs in the 
current solar system, (3) the annual flux of dynamically new lpcs inside 5 au, (4) the inner to outer oort cloud 
population ratio, (5) the derived origin probability ratio for 2000 0067, (6) the derived origin probability ratio for 2006 

SQ372- 



seen in this figure that as semimajor axes approach 10 3 
AU (the approximate value for SQ372) the Oort Cloud 
origin probability increases by ~2 orders of magnitude 
relative to the scattered disk This result supports the 
work of lEmel'vanenko et alj (|2005l ) who conclude that 
the Oort Cloud should efficiently produce a population 
of planet-crossing orbits with a > 500 AU beyond Sat- 
urn. 

Thus, Figure 2 indicates that 2006 SQ372 is 16 times 
more likely to originate from the Oort Cloud than from 
the scattered disk. We must keep in mind, however, that 
the results shown in Figure 2 depend on a number of 
quite uncertain parameters. Other combinations differ- 
ent from those used in Figure 2 are far from ruled out, 
and in Table 1 we list the resulting 2006 SQ372 origin 
probability ratios for a range of other possible normaliza- 
tion factors and inner Oort Cloud populations. It can be 
seen from this table that the probability of an Oort Cloud 
origin for 2006 SQ372 remains higher than the scattered 
disk for any combination and may be over 1000 times 
the probability of a scattered disk origin. From these re- 
sults we believe that 2006 SQ372 is the most distant LPC 
discovered to date. 

4.3. An Inner Oort Cloud Origin 

If 2006 SQ372 is from the Oort Cloud, even more in- 
triguing is that it almost certainly resided in the in- 
ner, unobserved regions of the Oort Cloud rather than 
the outer parts where all LPCs observed near Earth are 
thought to originate. The dynamical reasoning is as fol- 
lows. Planetary perturbations must have been responsi- 
ble for lowering the semimajor axis of this orbit to ^10 3 
AU, yet this random- walk proces s proceeds at a very slow 
pace. From lEhmcan et "ail (|1987l ). it can be shown that 
^ - r 12 / „3 XV2 

dt 
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where D is the average energy kick per perihelion pas- 
sage, do is the original semimajor axis, a/ is the final 
semimajor axis, and tjj is the time required for the semi- 
major axis to diffuse down to a/ from ao via planetary 
kicks. Setting a f to 800 AU and D(25 AU) to lO" 4 7 
AU" 1 (|Duncan et al.lll987t iFernandez fc Bruninill2000f l , 
we can then plot to as a function of initial semimajor 
axis in the Oort Cloud in Figure 3. As can be seen in 
this plot, tjj is 130 Myrs for a = 10 4 AU and increases 
for larger a . 




Fig. 3. — Plot of the seimajor axis and perihelion diffusion 
timescales vs. initial semimajor axis. The timescale for a peri- 
helion shift of 10 AU (the shift necessary to detach SQ372 from 
the planets) is shown by the dashed line, and the timescale for an 
orbital semimajor axis to random walk down to 800 AU is shown 
by the solid line. In addition, we have overplotted a histogram 
of the initial Oort Cloud semimajor axes for all objects somewhat 
similar SQ372 (20 AU < q < 30 AU, 10° < i < 30°, and a < 800 
AU) that we generated in our Oort Cloud simulations. 

Left unchecked, the planets would be able to draw 
down the semimajor axes of any Oort Cloud comet within 
a couple hundred million years. Before enough time can 
elapse for this to happen, however, most comets have 
their perihelia lifted back out of the planetary region 
again by the galactic tide. If a comet's perihelion is 
moved outward by —10 AU from 25 AU to 35 AU, the 
planets cease to have a large effect on the dynamics, and 
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the random-walk of semima jor axis stops. The t imescale 
for this to occur is given in Ibuncan et al.l (|1987f ) as 

, - , / Ag \ /25AU\ 1/2 /10 4 AU\ 2 

^=( L3Xl ° 7yr )(lOAljJ(— j {—) 

(4) 

Assuming a perihelion of 25 AU, t q is also plotted in 
Figure 3. To produce SQ372-like orbits from the Oort 
Cloud, to must be shorter than t q , otherwise a = 10 3 
AU will never be reached. As can be seen in this plot, 
this is only the case for a subset of Oort Cloud bodies. 
Only for initial semimajor axes below ^5000 AU can 
the semimajor axis be drawn down to 10 3 AU before 
the perihelion is lifted out of the planetary region. This 
process is verified when we examine the distribution of 
initial Oort Cloud semimajor axes for all of the SQ372- 
like orbits we generate, which is also shown in Figure 3. 
It can be seen in this figure that the distribution falls off 
quickly for a > 10 4 AU. Therefore, we can conclude from 
this argument that if 2006 SQ372 is from the Oort Cloud, 
it is very likely from the inner 10 4 AU of the Oort Cloud. 

4.4. Comparison with Other Objects 

4.4.1. Comet C/2003 A2 

Comet C/2003 A2 (Ideason et al.ll200l lGreenll200l is 
the only previously known LPC with perihelion beyond 
Saturn. Because its perihelion is still within 2 AU of 
Saturn (q — 11.4 AU), its dynamics differ significantly 
from 2006 SQ372. This object's proximity to the strong 
energy kicks of Saturn make it just as likely to come from 
the outer Oort Cloud as the inner Oort Cloud. This is 
because the timescale for Saturn to modify the semimajor 
axis of a comet is much shorter than Uranus or Neptune. 
As a result, our argument for the inner Oort Cloud origin 
of 2006 SQ372 cannot be applied to the orbit of Comet 
C/2003 A2, and we cannot constrain the region of the 
Oort Cloud where this comet previously resided. 

4.4.2. Sedna 

Given that we believe 2006 SQ372 is from the inner 
Oort Cloud, a comparison with Sedna, another p ur- 
ported inner Oort Cloud body (|Brown et al.l l2004T i. is 
justified. Sedna is thought to have been scattered on to 
its current orbit (q = 76 AU, a = 487 AU) by strong ex- 
ternal perturbations of the Sun's b irthplace environment 
early in the solar system's histor y ( Morbidelli fc Levisonl 
120041: iBrasser et alj|200d iKaib fc Quinnll2008fh Because 
external forces of this magnitude no longer act in the 
Sun's vicinity, Sedna's orbit essentially ceased evolving 
Gyrs ago. In contrast, because of a larger original semi- 
major axis, 2006 SQ372 has recently had its perihelion 
pushed closer to Earth by the Galactic tide and pass- 
ing stars. Thus, 2006 SQ372 represents a more proxi- 
mate (q > 15 AU) population of inner Oort Cloud bodies 
whose orbits are still continually evolving. 

4.4.3. 2000 00 67 

Out of all known TNO s, 2000 00 67 dMillis et all 
((2002h : iVeillet et all (pOOlh : q = 20.7 AU, a = 552 AU, 
i = 20°) has an orbit most similar to 2006 SQ372. When 
we use Figure 2 to determine the probability that it is 
from the Oort Cloud, we find that the Oort Cloud should 



produce orbits similar to 2000 OO67 at a rate that is at 
least 14 times that of the scattered disk, assuming the 
comet population numbers quoted in section 4.2. Fur- 
thermore, because its perihelion is far beyond the or- 
bit of Saturn, similar timescale arguments can be used 
to show that 2000 OOg7 must also come from the in- 
ner ~10 4 AU of the Port Cloud. Thus, as suggested in 
lEmel'vanenko et all (|2005l ). it appears that 2000 00 67 
is also a comet from the inner Oort Cloud. 

5. SUMMARY 

We have discovered an object (2006 SQ372) on an un- 
usual orbit with a perihelion of 24.2 AU and a semimajor 
axis of 796 AU. Because of the ability of the scattered 
disk and Oort Cloud to produce planet-crossing orbits 
in the outer Solar System we believe this object pre- 
viously resided in one of these regions. To determine 
which scenario is more probable we have simulated the 
production of similar orbits from both the scattered disk 
and the Oort Cloud. The results of our simulations in- 
dicate that similar objects are produced from the Oort 
Cloud anywhere from 2.3 to 1100 times more often than 
the scattered disk, depending on the the populations of 
both regions. We argue that a production rate ratio of 
16 or higher is most likely and that an Oort Cloud ori- 
gin for 2006 SQ372 is therefore most probable. Intrigu- 
ingly, an Oort Cloud origin for 2006 SQ372 implies that 
is very likely from the inner 10 4 AU of the Oort Cloud. 
Otherwise, its perihelion would have been removed from 
the planetary region before its semimajor axis could be 
drawn down to 10 3 AU. This line of reasoning can also be 
applied to show that 2000 0067 must have had a similar 
history. An inner Oort Cloud origin makes 2006 SQ372 
and 2000 OO67 unique with respect to all known LPCs. 
Furthermore, these bodies are members of an evolving, 
more proximate population of inner Oort Cloud bodies 
that are easier to detect because of their small perihelia 
compared to Sedna and other members of the extended 
scattered disk. 

Although we have shown 2006 SQ372 is a probable 
LPC, this object cannot provide much new information 
about the structure of the Oort Cloud by itself. The 
compilation of a large sample of similar Oort Cloud ob- 
jects passing through the outer planetary region can, 
however, place better constraints on both the size and 
structure of the Oort Cloud. The prospects for compil- 
ing such a sample in the near future are encouraging. 
2006 SQ372 was discovered by imaging the same fields of 
sky over multiple nights separated by periods of months 
during the SDSS-II SN survey. In terms of both limiting 
magnitude and sky coverage, this survey is modest com- 
pared to the large synoptic surveys such as Pan-STARRS 
ilKaiser et al.ll2002HJewittJl2003f ) and LSST (jlvezic et all 
120081 ) being planned for the coming decade. Thus, we 
can expect many objects similar to 2006 SQ372 to be 
discovered in the near future. Used in conjuction with 
dynamical modeling, these discoveries will provide many 
clues about the current size and structure of the Oort 
Cloud as well as the dynamical history of the solar sys- 
tem. 
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